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Abstract 

Dimensional reduction of finite temperature quantum field theories 
can be improved with help of continous renormalisation group steps. 
The method is applied to the integration of the lowest non-static (n = 
±1) modes of the finite temperature $ -theory. A second, physically 
important application is the integration of the Debye-screened Aq(x) 
static scalar potential in the gauged SU(2) Higgs model. 
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1 Introduction 



Thermal field theories at finite temperature can be equivalently formulated 
with help of countably infinite number of three-dimensional fields represent- 
ing each quantum field (Matsubara decomposition). In the three dimensional 
formulation each set of representatives can be classified at tree level as be- 
longing to the group of massive (m ~ T) or massless (static) components. 
(Fermions are represented exclusively by massive modes.) The reduction 
program consists of integrating out the massive components, their role being 
reduced in this way to influence (sometimes essentially) the interaction be- 
tween static modes. The fluctuations of these latter are widely accepted to 
govern temperature driven phase transitions in quantum field theory. 

The usual procedure of perturbative reduction [[l], ||] is of restricted valid- 
ity. Especially detailed quantitative non-universal characteristics (T c , order 
parameter discontinuity, surface tension, etc.) might be very sensitive to the 
correct dependence of the couplings of the reduced theory on the temper- 
ature and on the parameters of the original (T = 0) theory. For this one 
should achieve the most faithful mapping of the full finite temperature sys- 
tem onto the effective theory. In case of the QCD deconfining transition an 
answer has been attempted to this problem by matching numerically several 
observables of the full finite temperature theory to the values simulated with 
the three-dimensional effective theory || f|, [5|, ||. 

The use of an "exact" Renormalisation Group transformation for the local 
part of the action has been practised in critical phenomena already for two 
decades for the determination of the RG-flow and of the fixed point effective 
potential |7|, || [|, [K| [II], [12], [13|. To our knowledge, however, the investigation 
of its usefulness from the point of view of the dimensional reduction appears 
for the first time in the present paper. 

We shall test its applicability by completing the reduction in models where, 
after a partial (perturbatively performed) reduction only finite number of 
three-dimensional fields are left, but for some reason still a well-defined gap 
persists between their respective thermal masses. Intuitively, perturbative 
integration should hold for the most massive (highest) Matsubara modes, 
but considerable improvement is expected to result from a more accurate 
non-perturbative integration over modes, lying the closest to the lightest 
ones. 

The first question we should like to clarify is to what extent it is possible 
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to take into account higher loop effects in the reduction, which reflect the 
interaction between the non-static modes. Our proposal is to integrate out 
perturbatively only the modes with \n\ < n , and proceed to a second stage 
of integration with the couplings already modified in the first step. Actually, 
the application of the RG-approach in the second stage allows the realisation 
of a continous feedback of the modified couplings on the process of lowering 
the momentum scale for the fields to be integrated out. As a testing ground 
for the procedure outlined above, we shall investigate the finite temperature 
1-component $ 4 theory with the choice hq = 1. 

A similar situation is arrived at in the gauged SU(2) Higgs model, after 
having integrated out perturbatively all non-static modes. The phenomenon 
of Debye-screening creates a new distinct mass scale (m ~ gT), separat- 
ing the A (x) static modes from the rest |14|, [TI| |16j. It is then natural to 
continue with a second (hierarchical) integration of the 3 A -field compo- 
nents, and apply MC-simulation methods to the resulting three-dimensional 
gauge+Higgs effective system. Propositions for perturbative realisation of 
this program have been already published [|T7| , |T8| . 

In this paper we shall expose and test the application of the "exact" 
RG-transformation to situations corresponding to the two problems shortly 
outlined above. We are going to use a local version of this transformation, 
therefore effects of the field renormalisation will not be considered in the 
present paper. 

In Section 2 the general strategy of the application of Renormalisation 
Group motivated reduction will be discussed in detail. Thorough qualitative 
description of some intrinsic features of the procedure will be given, providing 
insight into the temperature range where the proposed RG aided reduction 
is expected to work. In Section 3 we introduce our test-models mentioned 
above (the lengthy, but interesting excercise of partial perturbative reduc- 
tion of the $ 4 -theory to the Uq = 1 level is relegated to the Appendix). In 
Section 4 the RG aided reduction will be realised on these systems under 
the assumption that the local potential density of the resulting effective the- 
ory can be well truncated to the linear combination of a finite number of 
lowest dimensional operators. The arising finite coupled set of differential 
infrared RG-flow equations will be carefully integrated, and the resulting 
couplings compared with the characteristics of the fully perturbative (some- 
times higher-loop) reductions. The summary of our results with a detailed 
discussion of the qualitative features of the final light theory is presented in 
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Section 5. 



2 The Strategy of the Renormalisation Group 
Aided Reduction 

The generic situation one deals with is a three-dimensional theory (arrived at 
in previous - possibly perturbative - reduction steps) with one light (0) and 
one heavy (<£>) field in interaction. At momentum scale A their interaction 
potential is known 

f/ A ($,0) = \/($,0,A,T) (1) 

and in particular one has an explicit expression for the corresponding mass- 
terms: 

|$= , =o = ™ 2 + HT 2 - b,AT ee Mj(A, T), 

0,^=0 = m " + a * T2 ~ & * AT = Ml(A,T). (2) 

These masses are renormalised from the four-dimensional (T = 0) point of 
view, but contain explicitly the induced three-dimensional "counterterm"- 
contribution depending linearly on the cut-off. If the reduction is realised 
with higher loop accuracy the linear cut-off dependence might be modulated 
by logarithmic terms. 

By assumption there is a relation between the values of the two bare 
masses at high momentum: 

Mj(A,T)<M|(A,T). (3) 

(m 2 is the mass parameter of the original theory renormalised by a condition 
imposed at T = 0.) 

The integration of the heavy field $ will be performed in infinitesimal 
steps, integrating in each step over its Fourier components lying in the thin 
layer k < |p| < k+Ak = K, Ak — > 0, and feeding back the resulting modified 
couplings into the next integration step. For an infinitesimally thin layer the 
1-loop contribution will be dominant. It has to be evaluated in a background 
of the fields 0l(x),$l(x) composed exclusively of Fourier components of 
momenta |p| < k. For the computation of the potential energy density 
(f)L, <3>£ can be set constant. 
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k^ = ~— fc 3 ln(fc 2 + ^($ L ,0 L )). (5) 



The matrix characterising the high spatial frequency fluctuations $#(x) 
at the Gaussian (1-loop) level is given by 

= Ml(K, T) + Cl (K, T)$| + c 2 (K, T)<f>l + & L <fc, <t>\) (4) 

The change in the potential energy is governed by the familiar differential 
equation 

dU k 1 
dk ' 4tt 2 

This equation should be integrated "downwards" on the momentum scale k, 
with the initial conditions (1),(2) fixed at k = A. 

For small field values of $^,0^ the quadratic form (4) is well approx- 
imated by the terms written out explicitly. Near the cut-off A the large 
negative contribution to M|(A,T), linear in A is compensated by adding in 
the argument of the logarithm to it k 2 = A 2 and the contribution to is 
real. In case M|(fc,T) would not run appropriately with k, the argument of 
the logarithm on the right hand side of (5) would change sign at a certain 
scale fc$(T), which would mean that a consistent integration of $ at that 
temperature would prove impossible. 

Since M$(k = 0, T) is the renormalised mass of the heavy field, for suf- 
ficiently high temperature one can expect with certainty it to be positive. 
Then the reduction can be performed meaningfully. There is a threshold 
temperature T$ at which M| crosses zero first as the temperature is being 
lowered. It presents a limit in temperature to the consistent reduction with 
the present simple formulation of the Renormalisation Group idea. It could 
well be that by appropriate fine-tuning of the cut-off procedure one could 
avoid the reduced potential becoming imaginary, in the same way as it has 
been demonstrated for the effective potential by ill , B| . However, this would 



make the actual procedure rather complicated and it seems to us not neces- 
sary for the purposes of the reduction on the basis of the following simple 
argument. 

We would like to study the immediate neighbourhood of the transition 
temperature T c . Therefore the reduction procedure is useful only if one finds 
T c > T$. It is plausible that the heavy-light mass-inequality (3) will persists 
all the way when rolling down on the momentum scale k: 

Ml(k,T) < Ml(k,T). (6) 
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Then repeating the above logical steps for the light degree of freedom we are 
led to the conclusion that the instability of the constant background field 
4>l, signalling the phase transition, will indeed occur at higher temperature, 
than T$. 

3 Finite temperature effective models with 
heavy- light field hierarchy 

The challenge we face is to find the simplest representation of the electroweak 
theory at finite temperature containing the smallest number of degrees of 
freedom, without loosing any essential aspect of the phase transition. The 
perturbative integration over n^O Matsubara modes seems to be intuitevely 
correct since their mass scale 2ttT is much larger than the scale characterising 
nonperturbative quantities, like magnetic screening. This latter is character- 



ized by the scale g T |20|, an d for small g it is even smaller than the "electric" 
Debye-screening mass-scale gT. 

The effective three-dimensional theory has been derived on the 1-loop 



level in different renormalisation schemes by different groups |fL8j , plf| . Fol- 



lowing [21] we use the form obtained with momentum cut-off regularisation 
in the scheme where the location and the value of the non-trivial minimum 
of the potential energy density of the effective model is kept at the values 
fixed on the classical level p2|. Its expression is the following: 



where 



S[A t , Ao, <P]=J d 3 x[L kin + U(A , <f>)], (7) 



L kin = + \[{di + igA^+m + igAi)<f>] + i(<Mo + ge^A\Al)\ 

fWA), 0) = \™W<\> + \m 2 D Al + §£ 2 A 2 0+0 + f(0 + 0) 2 

+ TiS(^o) 2 - [fU 2 , + (¥ + A)^ + 0]A, (8) 

-J = - 2 + (^ 2 + ^)T, rnl = lfT + ^t, (9) 
m 2 = m 2 (l - 3^((fs 2 + A) In + A In §§) - ^(45^ 2 + 20A + ^)), 
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A = A-^(^(|ln§-|ln0) + ^ln^ + 3(¥ + t) 2 ln^#) - 

-T^(¥ + ^ + f)- (10) 

Perturbative investigations of the phase transition have shown that near 
the transition point the screening mass of the A Q (~ gT) quanta is consid- 
erably larger than the masses of the Higgs and of the magnetic gauge field 
fluctuations. For this reason we assume (and shall test subsequently the con- 
sistency of the assumption) that in this model the A multiplet represents 
the heavy fields which can be integrated out in some careful procedure. 

An analogous situation can be created by construction also within the 
familiar 1-component $ 4 -model, if instead of integrating over all non-static 
fields at once, as a first step only the \n\ > 2 modes are integrated out. As 
a result an effective theory with a light (M ~ AT) and a heavy (M ~ 2nT) 
field arises. (One should note the difference in the size of the heavy and light 
mass-scales between the two models). 

We expect the subsequent integration of the heavy field will display also 
some features of the 2-loop reduction, but not any really dramatic deviation 
from the standard 1-step reduction. Applying a non-perturbative integration 
technique to the heavy component one can test the stability of the standard 
reduction (in the weak coupling regime) against this expectation. 

The derivation of the effective heavy-light theory in the <3> 4 -theory repre- 
sents a technically quite involved perturbative procedure. Since in the main 
text we wish to concentrate on features of the non-perturbative reduction, 
we describe details of the perturbative part of the reduction in Appendix 
A. The result of this (partial) reduction is given by the following effective 
action: 

Ss D [(f>o, 01] =/^{|(W + \W + Wl + 2|0i| 2 )K + if: - 0) 

+ (27rT) 2 |0 1 | 2 + A[04 + 12(1 + aA )02| 01 |2 + 6(1 + 6A) | 0i |4] 

+ Wo,0l]}, (11) 

where 

00 = ^0, </>i = yJp*i, A = AT, (12) 
and the numerical constants (see Appendix A) 

a = 3.56xl(T 3 , b = 3.17 x 1(T 3 . (13) 
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The normalisation condition required at this level to ensure the 4-dimensional 
ultraviolet finiteness of the theory is of the following form: 



^rU=o-A. (14) 

This condition explains the presence of finite (C(A 2 )) corrections to the terms 
|0! | 4 and |0i| 2 0q above. 

Also the derivation of the sixth-power contribution to the potential is 
given in Appendix A. It is of the form: 

AU 6 = AU^ + AUf ] (15) 

with the following explicit expressions for the terms on the right hand side: 

AUP = -^C^l^l^ + ^| 6 ^), 

AC/ 6 (2) = 5.40183 1(T 6 M 6 + 2.7566 l(r 5 A 2 0g|0i| 2 M 2 

+6.6586 1O- 6 A 2 |0i| 4 M 2 - 1.6991 lO^A^ 2 ,^ 4 (16) 

(M 2 = m 2 + A(0o + 2|0i| 2 )/2, u n = 2imT). 

This contribution is usually neglected, since it is suppressed at high tem- 
perature by $ 2 /T 2 relative to the quartic terms, when |$| «T. It might 
prove interesting to include it into the non-perturbative treatment of the 
0i-integration, at least to check the validity of the truncation at the quartic 
term. 

In the next section the evolution equations for the potential arising from 
the non-perturbative integration over subsequent momentum layers of the 
0i field will be worked out explicitly. In the approximation, when the 0o- 
potential of the resulting theory is projected onto a quartic polynomial, 
these equations can be exactly mapped onto the flow equations of the Higgs- 
potential of the electroweak theory induced by the ^-integration. This al- 
lows a uniform formal treatment of the two physically different heavy-light 
systems from the point of view of the reduction. 



4 Non-perturbative integration of heavy fields 
(finite polynomial approximation) 

The detailed analysis of the non-perturbative integration will be illustrated 
with the 0o _ 0i system. At the end of the section we describe the corre- 
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sponding results for the effective theory of the electroweak phase transition. 

The first step of the procedure consists of the application of two subse- 
quent Hubbard-Stratonovich transformations, through which the operators 
~ |0i | 6 and ~ |0i | 4 can be expressed in a form apparently quadratic in 0i 
with help of two auxiliary Gaussian variables {xi,X2)- Next, the integra- 
tion in an infinitesimal momentum layer over the high spatial frequency part 
of the 0i-field is performed in the background of the slowly varying fields 
0ol, 0il, Xil, X2L- The accuracy of the result can be optimised by choosing 
appropriate 0o, 0i dependence for the auxiliary fields xi an d Xi- 

The result of the integration is of non-polynomial nature. Due to the 
large mass of the n = 1 mode the interaction is short ranged and screens 
the infrared behaviour. In other words, the scale dependence is given by 
the ultraviolet scaling relations for all length scales. This allows to restrict 
the solution for up to the sixth order terms of the local potential and to 
expand the renormalization group equation in the interaction, (p8|). By this 
argument throughout this paper we approximate the full potential energy 
density by a finite order polynomial (namely, up to 4 or 6 ). The result of 
an infinitesimal step can be expressed then in terms of a set of coupled first 
order non-linear differential equations. They describe the flow of the theory 
in a restricted coupling space as the momentum scale bounding the support 
of the heavy field (0i) is lowered. 

The effective light (0o) theory is arrived at when the above bound reaches 
zero, that is the 0i integration being completed. It is worthwhile to empha- 
size that 0i does not fully vanish from the theory. It is present as a constant 
background, the integration of the 0o-theory (for instance the calculation of 
its effective potential) should be performed in such background. This proce- 
dure clearly will react back on the final 0i dependence of this two-variable 
function. It might look just a curiosity, that our two-step integration over 
the non-static modes allows us searching for a non-trivial minimum in a two- 
variable (0o, 0i ) domain. Usually it is correct to assume that the absolute 
minimum lies for any temperature in the 0i = slice. Still, it could occur 
in some segment of the coupling space, that a ground state which breaks 
invariance under r-translation is realised. 
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4.1 The generalised Hubbard-Stratonovich transforma- 
tion for the $ 4 -theory 

In the following analysis the projection of the potential energy density will 
be restricted to dimension-6 operators. The starting action at momentum 
scale K is parametrized in the form: 

£[00, = S^[\ipi<h? + \d*fa\ 2 + iM K )^ + iiMKWo 
+((2ttT) 2 + B^K))^ 2 + ±B 2 {K)\^ + C{K)<ft\k\* 

+D 1 (K)$ + D 2 {K)4>i\<t>i\ 2 + DsWffilfal* + D A {K)\<j> x \% (17) 

A two-step generalisation of the Hubbard-Stratonovich transformation will 
be used, based on the following simple identities: 

/££ d xi eM\x\ - £iXiW 2 - 6xi|0il 4 + \Q\4n\ 4 ] 

~exp(-66l0i| 6 + 0(l0i| 8 ), 
/St^ 2 exp[ix!-e3X2|</'i| 2 ] ~exp(-|ell0i| 4 )- (18) 

Requiring the relations 

66 = £4, ^ 3 2 = + ^ + 6xi - (19) 

one can use the equalities flTBp for rewriting the the action in extended form 
(possible terms of O(|0i| 8 ) are omitted): 

S[MiU = J d 3 x{±{d i( f>o) 2 + l^xl 2 + |A!(K)02 + iA 2 (if)^ 

[(2nT) 2 + B 1 (K) + £1X1 + 6X2 + C(lf)0g + D{K)<^\\<k\ 2 •} (20) 

For the ^-integrals we follow the usual approximation of saddle point 
dominance. The long wavelength parts of both 0o and 4>i are replaced by 
constants, which enables one to perform the integration over the short wave- 
length part of 0i easily: 

^ jf* dp p 2 \og(p 2 + (2ixT) 2 + B x + xi + X2 + C4> 2 + D 2 <f>i). (21) 
For the scaled quantities 

Xi = 6Xi, £2 = 6X2 (22) 
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one gets the following gap equations: 

Xi = (£! 2 + iW£!)(l0ii 2 + 'i(A)), x 2 = £ 3 2 (l0i| 2 + A(A)). (23) 
Here the notation A = xi + X2 + B\ + C<t>o + -^2^0 an d the integrals 

1(A)-— l K t*L f 2 4) 

n{ ' 2tt 2 A (p 2 + (2vrT) 2 + A)" 1 ' 

have been introduced. 

In the classical extended action (pOD, evaluated with the long wavelength 
parts of the fields the saddle point contribution can be shown to differ from 
the classical expression -^3^0 1 ^1 1 4 + -Si|0i| 4 /12 + -D/tl^iJ 6 only in terms pro- 
portional to higher powers of I\. They give negligible contribution to the 
couplings at scale k, when the interval (k, K) becomes infinitesimal. 

Therefore, the running of the couplings comes exclusively from (|2l|) . Its 
field-dependent part can be expanded into powers of |</>i| 2 and 0q : 

AU = J2 -^—(Xi + X2 + C<t>l + D 2 ^) n / n (S0. (25) 

n=l U 

Since each term is multiplied by an integral over the infinitesimal range 
(k,K), for the differential running that part of Xi + X2 which is propor- 
tional to another similar integral I n does not play any role. Combining the 
two equations of ( [2"3"D and the definitions of £ 2 , £f, one can show that 

Xi + X2 = 2 J D 3 0ol0i| 2 + 3AI04I 4 + ^ 2 |0i| 4 + Od^l 6 ) (26) 

o 

is the relevant (classical) value to be used for the derivation of the differen- 
tial running equations of the various couplings. The (9(|0i| 6 ) contribution 
would produce a linearly divergent term in the running of -D4. It should be 
omitted by the neglect of the irrelevant |0i| 8 term. The full set of evolution 
equations appears in Appendix B. Here we discuss analytically the evolution 
of dimension 2 and 4 operators. 

When one projects the potential onto a fourth order polynomial, it is 
sufficient to keep only the first two terms of fl2"5|). The most convenient is to 
introduce dimensionless variables by scaling all three-dimensional couplings 
by appropriate powers of 2tiT: 

k „ _ Ai u B\ 



2wT , "1 — ( 2 7rT) 2 ' U1 ~ (2tvT) 2 ' 
2-kTi U2 2ttT' u 2-kT- 



a-2 - 5=r, b 2 - c - =g*. (27) 
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Then the following simple set of coupled first order differential equations 
arises: 



dai c x 2 dbi b 2 x 2 

dx ~ k 2 x 2 +l+bi ' dx ~ 12-n 2 x 2 +l+b 1 ' 

dc &2C x 2 da,2 6c 2 ™- 2 



dx 12-k 2 (x 2 +l+b 1 ) 2 ' dx TT 2 (x 2 +l+b 1 ) 2 ' 

(28) 



db2 _ 2:2 
"dx ~~ 12k 2 (x 2 +l+b 1 ) 2 



Comparing the right hand sides one easily finds some natural relationships 
between the evolution of various couplings: 

c = Qih, a 2 = nQ% + Q 2 , ai = 12Qi&i + Q 3 , (29) 

where Qi are integration constants. These relations imply that the true task 
is the solution of the coupled evolution equations of b\ and b 2 . 

In the weak coupling regime, where also b% can be assumed in pertur- 
bative sense to be small, in the first approximation b% can be neglected in 
the denominator of the right hand side of the equation for 6 2 and can be 
represented by an expansion of the denominator up to linear terms on the 
right hand side of the equation of b\ (One works to quadratic order on the 
right hand sides of (p8|)). Then, it is easy to find explicit analytic solutions: 



b 2 (x) 



24tt 2 (1+x 2 ) 



x+(l+x 2 )(247r 2 Q 4 -tan- 1 x) ' 
= , + (l + , 2 )(24:4-tan-^) ^5 + 2^ *-*)]. (30) 

The integration constants can be found from the intial conditions which are 
the couplings at momentum scale A, arising from the 1-loop integration over 
the |n| > 2 modes: 

ai(x A ) = bx(x A ) = { ™ T)2 + ^2 - 
a 2 (x A ) = ^, b 2 (x A ) = £(l + b\), c (x A ) = A(i + aA ). (31) 

The constants Qi are found with straightforward algebra when equating (|3T|) 
to the large x A limiting values of the corresponding functions in (p9|) and 

The couplings reflecting the effect of 0i-fluctuations are given by ai(0), 6j(0), 
c(0). The couplings a i: c determine the effective fluctuation dynamics of the 



12 



static modes (0o)- h characterize the effective potential of constant <\>\ config- 
urations. The assumption of the smallness of B2 and C which was essential 
for the consistent truncation of the infinite set of RG-equations to ( pgp is 
ensured by restricting A to small values. Then, it is sufficient to give the first 
few terms of the power series of different couplings with respect to A: 

6a(0) « §(1 + A(6 - 3^)), c(0) « £(1 + X(a - ^)), 
«a(0)«£(l-jg*), 
6i(0) « ^(1 - 3^(1 - A(^ - 6))) + + ^ + 6) 

«l(0) ~ ( 2 ^T) 2 ~~ T^( 1 + _ 3^))) + 96^ _ &?(l92^ + 8^ + fl ) 

+ 8^(-| + A («+T^))- (32) 



4.2 Justification of the truncation 

The solution of the approximate renormalization group equations ( p8|) is 
based on the assumption 1 + x 2 > \b\\ which amounts to the inequalities 

m 2 13A 3AA A 2 , , 

+ 7^ - TcZAm < 1 + 7^wTn2' (33) 



(2ttT) 2 96vr 2 16vr 4 T (2vrT)< 

or 



m 2 13A 3AA 1 A 2 

> - 1 " T^vr (34) 



(2vrT) 2 96vr 2 16tt 4 T (2ttT) 2 
(|33| ) can be taken for granted since its violation requires 

°< r2 (^-f -|£)<» 2 . (») 

and A = O(AT). 

(|34T) represents a more stringent condition, 



m » > _W^-^-i^V (36) 
V 24 56ttV v ; 

For m 2 = — 0((27rT) 2 ) the <\>\ system becomes critical for O = 0. But this 
infrared instability is misleading, it comes from the parametrization of the 
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effective theory only. The coupling constants ax, a>2 and c were obtained 
by expanding the effective potential around 0o = 0, so they are relevant in 
describing the fluctuations in the symmetrical phase or in the vicinity of the 
phase transition. For deeply broken symmetry systems the mass term C0q 
stabilizes the 0i system. In this case the usefulness of the integration of 4>i 
alone is questionable and the O ~ <j>i systems should be treated together. 

These arguments demonstrate that the numerical solution of the original 
equation ( p8"D would not provide in the weak coupling regime quantitatively 
different information on the reduced action compared to the approximation 
where the denominator is expanded in powers of b±. 



4.3 Comparison with results of simplified treatments 
of the (^-integration 



The coefficients a, appearing in fl32|) can be compared with less sophisticated 
ways of integrating out 4>i- 

The simplest is the 1-loop integration with constant (fro background, which 
is a version of the hierarchical integration advocated previously within the 
3D effective theory of the electroweak phase transition ||17|| . The quadratic 
form in the Fourier-space is simply determined by 

k 2 + m 2 + (2nT) 2 + + M (i + aX ), (37) 

which leads to the following contribution to the effective 0o _ theory: 

^(1 + aX) 1 -^ - £(m 2 + (2vrT) 2 + ^ + M(i + aX))W (38) 

For A0Q «T 2 one expands this expression in power series with respect to 
A up to terms C(A 2 ), what is equivalent to retain terms only up to 0{4>q). 
In the expressions below also the leading 0(m 2 /T 2 ) corrections are included. 
The results for the coefficients the following: 

Ol(0) = (2^)2 (1 - i^2-(l + a X) + 30727r 4 ) + ~ + lsi!?) 

a 2 (0) = A(i__^ A (l-^)). (39) 
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One instantly recognizes that the difference between (|32| ) and ( |39| ) of the 
resulting couplings starts at 0(\ 2 ). 

One might improve upon the simplest 1-loop hierarchical integration by 
applying first the HS-transformation, as described above. Essentially this 
amounts to the evaluation of /„ in (|24"D with k = 0, K = A, and keeping 
only the <fto background different from zero. The result is 

a l(0) = ( 2 ^r)2 (1 - 1^(1 + + 3072tt 4 ) + 9^2 ~~ ^2"( a + 19^2 ) 

-8^(|-A(a-3^)), (40) 

while the expression of 02(0) agrees with the 1-loop (and also the RG) result 
toO(A 2 ). 

The change in the shape of the effective $0 — $1 potential depends on 
the renormalisation conditions. Since they were imposed at an intermediate 
stage of the integration, it is difficult to compare the above result with, for 
instance, an 0(X 2 ) evaluation of the effective potential without reduction. At 
present the best one can do is to evaluate the phase transition temperature 
a la Landau from the vanishing of the finite part of ai(0) in the different 
approximations. The original estimate of Linde and Kirzhnitz receives 
an O(X ) correction as follows: 

A -T 2 C = 1 + 0.0596A (RG), 1 + 0.1187A (pert), 1 + 0.0427 \(s addle), 



24m 2 ' 

(41) 

as one calculates it from eqs. (|32| , |39| , pOD , respectively. This can be compared 



with the critical temperature calculated in the same way from the effective 



static $0 theory arrived at with the reduction peformed on 2-loop level [^4 
1 + 0.0462A. This shows that the saddle point integration in the present 
context is "equivalent" to the 2-loop calculation, while the non-perturbative 
integration "sums up" the two- and higher loop contributions into a reason- 
ably nearby estimate. 

4.4 Non-perturbative Aq- integration in the effective elec- 
troweak theory 

The integration of the non-static modes has provided a new distinct screening 
scale, the Debye-mass of A , (pf). This mass has been shown in perturbative 
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treatments to be consistently larger, near the phase transition point, than the 
effective masses of the magnetic vector- and Higgs-field fluctuations. Based 
on this observation a treatment fully analogous to that of $1 can be envisaged 
also for A . 

For the simplicity of presentation we restrict the potential energy density 
to its quartic projection. In Appendix B we show, that the results are stable, 
when the projection is extended to dimension 6 operators. 

Then the Lagrangian of the theory at scale K can be parametrized as 



L = \FijFij + |(A0) f (A0) + §(AA)) 2 + \A x (K)tf<l> + ±A 2 (K)(tf<i>) 2 
+ \B X {K)AI + \B 2 {K){Alf + C{K)^ct>Al (42) 



After the HS-transformation is applied to the quartic A -potential, the La- 
grangian density gets the form: 



The integration over the Fourier components of A belonging to the (k, K)- 
layer contributes to the potential energy density 



The omission of the magnetic vector potential from the background, in which 
this expression has been evaluated corresponds to neglecting wave function 
renormalisation effects. 

The optimisation of the HS-transformation is achieved through the choice 
of x according to the gap-equation: 



It has the following iterative solution to quadratic order in the couplings: 



Xo = 2B 2 (A* + 3/0, xi = -4£ 2 / 2 (2C0V + 2B 2 {A 2 + 3/0), (46) 



/ = \FijFij + |(A0) f (A0) + i(AA)) 2 + \A^H + 

-^- 2 X 2 + \{B l + x)Al + C^^Al 



(43) 




(44) 




(45) 



where the notation 




(47) 
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has been used. 

In this case it seems to be more convenient to scale dimensionfull quan- 



tities by appropriate powers of m D = y5/6gT. On this scale we expect 
b\ = Bi/m 2 D — 1 + O(g). With this assumption we break up 61 into 1 + b\. 
bi is expected to sum up higher order contributions, but not to exceed the 
order of magnitude of g. Repeating the same steps as in case of the 4 theory 
one arrives at the following set of renormalisation group equations: 

dai 3c x 2 db± 3fe2 x 2 



dx it 2 x 2 +l+b± ' dx ir 2 x 2 +l+bi ' 

da 2 _ 36c 2 x 2 db 2 _ 3b 2 . x 2 

dx -it 2 (x 2 +l+bi) 2 ' dx it 2 (x 2 +l+bi) 2 ' 

dc 3&2C x 2 I AQ\ 

dx 7T 2 (x 2 +l+b 1 ) 2 - 

Like in the previous model, again some natural relations can be found be- 
tween the quadratic couplings, as well as the quartic ones: 

c = Q 1 b 2 , a x = Qxfei + Q 2 , a 2 = 12Q% + Q 3 , (49) 

therefore the problem is again reduced to the solution of the coupled set of 
the two equations for bi and b 2 . With the replacement 

62 - 1 (50) 

in (|48| ) one finds a form of the differential equations for 6 1; b 2 which is identical 
to fl28l) . It has the same solution with some constants of integrations, which 
can be determined from the initial conditions set at x A = A/m^: 



Oi(^a) - 1(7^2 + (4 + T^a) - J! (f + ^)f^")' M^a) - f^T 4 



~{gT) 2 1 V16 1 I2g 2 ' V 6^4 1 g 2 > 2tt 2 h "XV" AJ " 2n 2 

a 2 (x A ) = x fll, c (x A ) = 4§, b 2 (x A ) = Jl^. (51) 



The analytic solutions are given with help of (p0|), if on the basis of the 
assumption b x « 1 the j4 -propagators in ([48]) are expanded in powers of bi, 
keeping on the right hand sides terms up to quadratic order in the couplings. 

After some ugly algebra one arrives at the following expressions for the 
effective couplings arising from the RG-improved integration over A : 



62(0) - b 2 (x A ) 



289g b 1 

1024?r 5 1 , fa 21s. 



H 

" a i287r J 
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c(0) - c(x A ) = 



51.9 4 L 



2560vr J i I ft 173^ 
1+ V 5l28^ 



a 2 (0) - a 2 (x A ) = -gg - ^ 17g 3 - 



Mo) - B^) = v/f + ^|g)x A - f(i + ^g^)]—^ 



5 128ir^ 



ai (0) - ai (x A ) = + J\ g )s A - f (1 + Jfg)] - * , (52) 



5 128ir a 

In all couplings in the physical g-region (~ 2/3) the common denominator 
can be omitted since it gives ~ 1/1000 relative correction. 

The expression of b\ receives 0(g 3 ) corrections, justifying the starting 
assumption of b\ « 1. In the common square bracket of the expressions 
of ai(0) — ai(x\) and &i(0) — bi(x\) the O(g) corrections are of cca. 15% 
relative importance in the interesting g-range. 

For comparison we quote the corrections to the potential energy density 



from the one-step 1-loop integration of A JTy j : 



( fll (0) - ai (x A ))[l - loop] = -^1^ + VI^XA, 

MO) - a 2 (x A ))[l - loop] = (53) 

One recognizes that the leading corrections to the scaled mass parameter, as 
calculated with both methods, coincide. The non-trivial result is the next 
term, which appears only in fl52|) and mimicks the effect of a reduction on 
the 2-loop level. We conclude that the improved integration leads to higher 
order differences in the couplings. They might influence the result of the 
numerical simulation of the resulting gauged Higgs system appreciably, since 
their contribution represents a 10-15% correction relative to the 1-loop A - 
integration. 



5 Discussion and conclusions 

The derivation of an effective theory for three dimensional static fields is 
presented in this paper by means of different approximations. The idea of 
using effective models originates from combining different methods in dealing 
with different degrees of freedom. 

At high temperature the natural strategy is to eliminate first the non- 
static modes with the simplest one-loop approximation. The more delicate 
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problem of infrared sensitive static modes is dealt with in the effective model 
where the dynamics of the non-static modes is incorporated into the effective 
vertices. Even the simplest leading order perturbative solution of the effective 
theory represents a resummation of the thermal mass. 

This well-known scheme is a rudimentary application of the scheme of 
the renormalization group where the contributions of the modes which have 
been eliminated are taken into account in the next elimination step. In 
order to improve the effective model we pursued more systematically the 
renormalisation group strategy in the elimination of some massive non-static 
modes in the 1-component scalar model and of the Aq components of the 
vector fields in the SU(2) Higgs model. 

In the scalar model the n > 1 Matsubara modes were integrated out in the 
first step in the one-loop, independent mode approximation resulting in an 
effective theory for the and 0o modes. By this we ensure that in the second 
step when the n = 1 modes are eliminated the interactions between the n > 1 
and the n — 1 modes are retained. Since 0i is the lightest heavy mode 
we performed a partial resummation of the perturbation expansion during 
its elimination. This consisted of applying a saddle point approximation 
for a Hubbard- Stratonovich auxiliary field representing nonlinear it ies of the 
ineratctions of the 4>\ field, and solving the renormalization group equations 
for the relevant operators. 

Higher loop contributions in the solution (|30|) are generated by irrelevant, 
higher order operators when projected back onto the relevant ones during 
the one-loop elimination step. For the solution of the renormalization group 
equation such UV divergence structure is expected which is different from 
that of a systematic multi-loop elimination. The difference between the one- 
loop and the renormalization group solution manifests itself in the higher 
loop contribution, > i n (p2[) - 

In the framework of the renormalized perturbation expansion all UV di- 
vergent contributions are treated as small perturbations because they are 
multiplied by a positive power of the coupling constant. Our renormaliza- 
tion group refers to the bare theory and does not distinguish between the 
renormalized part of the lagrangian and the counterterms. Thus it is not 
obvious to what extent the results obtained with the help of this equation 
would correspond to a partial resummation of renormalized graphs. 

It turns out that the correspondence between the partial resummation of 
the renormalized perturbation expansion and the solution of the renormal- 
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ization group equation can actually be maintained. To see this we start with 
the observation that the scheme of the renormalized perturbation expansion 
could be violated if the inequality 

> \m 2 + 13AT 2 /24| (54) 

is fulfilled, i.e. when the mass counterterm is larger than the finite part. For 
m 2 < and A ~ this gives 

x- — >-— (55) 

Since x 2 > x > 1 one can treat the counterterm as small beside the kinetic 
energy and the expansion in the counterterms and the finite parts is justified 
by the same stroke. (This result holds in four dimensions, too, since the mass 
counterterm, though being quadratically divergent, is suppressed by A when 
compared with A 2 .) 

According to the argument presented above the UV divergent parts can 
be treated perturbatively for A << 1. Since A has only finite renormaliza- 
tion this inequality remains valid in the UV regime. Still, the "small" UV 
divergent part of the solution of the renormalization group equation makes 
the renormalization of the theory more involved. If the effective theory is 
derived and solved in the n-loop approximation then the general proof of 
the renormalizabiltiy allows to eliminate the cut-off. When the derivation 
and the solution of the effective theory are based on different approxima- 
tions then the divergent parts of the bare coupling constants of the effective 
theory do not cancel the divergences in the solution. In this manner unless 
the derivation and the solution of the effective theory match the resulting 
model is, strictly speaking, non-renormalizable. 

This can be illustrated by the inspection of fl32|) and (0). The one-loop 
integration of <\>\ generates ultraviolet divergences in (R8|) which are exactly 
cancelled by subsequent one-loop integration of (fi . But when, inconsistently, 
tree-level approximation is used for 0o then the effective potential remains 
cut-off dependent. In an analogous fashion the "bare" mass, ai(0) for 0o in 
(p2|) was obtained by the renormalisation group improved perturbation ex- 
pansion and its cut-off dependence will not be compensated for by its straight 
or any improved one-loop integration. In the same time the mass of 0i, &i(0), 
displays different ultraviolet divergences which remain unchanged after the 
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integration over 0o- At the end one finds different cut-off dependence for the 
masses of 0o, 0i and these infrared quantities can not be kept simultaneously 
finite. 

A pragmatic way out from this problem is offered by the careful exami- 
nation of the cut-off dependences. If the derivation and the solution differ in 
higher order of the perturbation expansion then the mismatch is 0(X n AT), 
n > 1. The cut-off should be large enough compared to the light mass scale, 
mi = J m? + - 13 2 J 2 , and small enough to keep this mismatch smaller than 
the characteristic heavy mass square, rrih = 2nT. If the expected scaling 
behaviour can be established for 

rm < A < (56) 

then the given solution of the effective theory is acceptable in the scaling 
window (jj6|). 

It might happen that in the relevant cut-off range only a somewhat more 
modest approach can be followed, corresponding to a phenomenological scal- 
ing law for the mass squares, 

m 2 {k = 0,T,A) = m 2 {T)+TAf{T/A). (57) 

If such relation is found to be fulfillled over an extended cut-off range then 
one might conjecture that it belongs to a certain renormalizable solution of 
the full theory. 

This incertitude in the use of the divergent parts of the effective theory is 
also a major problem in extracting physical results for the finite temperature 
phase transition of the SU(2) Higgs model when its 3 dimensional effective 
representation is being used, for instance, in the calculation of T c . The 
unsealed form determined in this paper looks like 

Seff = + I(A0) f (A0)| + m2(T)0t + JjA 3 (T)(W], 

ml(T) = m *+ [±g> + ±X- fgyf (1 + v/faO— ^]T 2 

AT [9„2 , \ _ 3 2/i I /6 5gN 1 1 

2^[j9 j9 U + V 5 8^J 1 7f j7g- 



+ V 5l5i^ 



HT) = [X-J I M)-^]T. (58) 



i+ 



5 1287T 3 
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It is very important for non-perturbative integrations of this effective 
model to test, how well the expected cut-off dependence (displayed in the 
second line of the expression of ml) fits the actually observed scaling be- 
haviour. The theoretical inaccuracy of the result is assessed by discussing 
the sensitivity of the critical data (T c , etc.) to the 0(g A T 2 ) and 0(g 3 AT) 
terms in m 2 (T). 

Appendix A 

1-loop integration of the \n\ > 2 Matsubara modes 

The 1-loop integration over the n^0,±l Matsubara modes with space- 
independent, r-dependent background 



starts by writing explicitely the action up to terms quadratic in the \n\ > 2 
Matsubara modes: 



S = (3V[\m 2 {<S> 2 + 2|$ 1 | 2 ) + £A($g + 12$g|$i| 2 + 6|$!| 4 )] 
+/3[X + f($ 2 + 2|$ 1 | 2 )]E M 0^0M 

+P^0 E' fc ,„($!0fc,„0fc,n+1 + $10 fe ,„0* fcin+1 ) 
+ P\ E M ($?0M0M + 2 + <t>l,n<t>k,n + 2) 
+ ^E[3$ o $?0O,-2 + 3$ $f 00,2 + $?0O,-3 + $f 0O, 3 ] + W 3 )(A2) 



From the primed sums the modes n — 0, ±1 are always omitted. 

The matrix of the quadratic form splits up into two identical disjoint 
blocs for fixed k, one for Matsubara modes n > 2 and the other for n < —2: 



r) = $(r) + <f>(x,r) 
$(r) = $ + $1 exp(icur) + $^ exp(-io;r), 

r) = En^o,±i 4>n{x) exp(iumr). 



= 2ttT, 



(Al) 



/ A; 2 + cu 2 + M 2 
A$i$ 



A$*$ 
A; 2 + luI + M 2 

A$i$ 








...\ 



M(k) 



A$2 

2 1 

V o 



A<I>i$o 

k 2 + to 2 + M 2 



... 
... 
/ 



(A3) 



where M 2 




($ 2 + 2|$!| 2 ). 
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The contribution to the potential part of the effective action is of the 
form 

£ log det M(k) - /3N* n M(0)^N m , (A4) 



with 



N n = \VV 



Id>*3 

5 1 



V 



(A5) 



The argument of M in the second term of ( Pjj ) tells that the inverse of the 
matrix ( |A3| ) should be evaluated at zero spatial momentum. Since we are 
interested in the reduced action up to terms of C , ($ 6 ), it is sufficient to keep 
only the diagonal terms in this inverse matrix 



M(O)- 1 



/ (m 2 + ujfj 



V 



2\-l 







(m 2 + lu 2 ) 



2\-l 



(A6) 



The contribution of this term to the sixth power part of the potential is 
calculated readily: 



(i> 



+ 



m 2 + u>2 36 



l$i| 



m 2 + uj 1 



(A7) 



The expression of the first term in 
terms in $o and $i reads 



expanded up to sixth power 



£ fc logdetM(^E fe , n {^ 
-A 2 ^! 2 l - ^ 4 



M 4 



M 6 



i 2(fc 2 +C 2 ) 2 + 



3(fc 2 +aj 2 ) 3 
A 3 ^ 2 |fi| 4 



+A 2 $ 2 1 |$i| 2 M 2 ' 
+^|$ 1 | 4 M 2 ( 



i 



+ 



i 



.(fe2+o;2)(fc2+^ +1 )2 "T" (fc2+o;2)2 (fc 2 +w j. + i) ; 

r)}- 



+ 



(k*+u>l)*(k*+u,* +2 ) ^ (fc 2 +^)(fc2+^ +2 ) 

Direct evaluation of the contributions proportional to $ 2 leads to 
1 A 2 1ST 2 SAT 

a^^ aK + 2 | 4i | 2)( -^ + HL 3 ^ 



4tt 2 



(A8) 



(A9) 
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For the evaluation of the fourth and sixth power contributions from 
it is very convenient to make use of the "mixed" or "Saclay" -representation 
of the n-sums P51: 



m = ^J^fPdT(-M± _ Ad^lV- - ^\^e 2 ^)G'{r) 2 
Af/ 6 (2) = (3V J £k Jf dnG'in) fg dr 2 G ' {r x + r 2 )C (r 2 ) 



~3 U 



x[^f + A 2 $ 2 |$i| 2 MV" Tl + ^ 1 \ 4 M 2 e 2iujTl 



- A 3 $2|$ 1 |V aj ( Tl+2r2 )]. (A10) 
Here the mixed finite temperature propagator is of the form 

G\t) = T£ j¥0)±1 e-^^ = G(r, k) - * - p^(e^ + e— ) 

G(r, fc) = i E s =± / s (fc)e"^, < r < /?, 

/- = n fc> / + = l+n fc , n fc = (e" fc - I)" 1 . (All) 

We elaborate on the computation of the quartic part, where the following 
relation, valid for any / 7^ can be exploited: 

fP J^piluTnif-Yi _ l+2n fc 2T I 1 , 1 N 

j am lt yi ) — k ^ k 2 + p^ fc2+a;2 V k 2 +(l-l) 2 u 2 ^ k 2 +(l+l) 2 uj 2 > 

~ k*(J+lW) + fc*(Jw)( 5 M + ^.-l) + (k 2 +u> 2 ) 2 (^.2 + 5/ _ 2 )- (A12) 

Upon the application of ( A.12 ) in the expression of AC/4 one arrives at the 
following representation containing only 1-variable integrals: 



AC/ 4 = -^(^ + 2|$ 1 | 2 ) 2 /^[^(l + 2-^) + 



2iE 2 (e a: -l) 2 



(x 2 +47r 2 ) 2 ^ 

A 2 



-(IV^^ 2 j dxx 2 [^^(l + ^) - , 2( , 2 2 +47r2) 



(a; 2 +47r 2 )(x 2 +167r 2 )J 

~~ |$l| 4 / [4-^2-^5) (1 + 2^zy) — a;2 ( x2+167r2 ) — (x 2 +47r 2 ) 2 

~~ (x 2 +47r 2 )(a; 2 +367r 2 )]- (A13) 

Here the notation x = (3k has been introduced, and for the divergent integrals 
a sharp cut-off A/? is understood. 
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From the large x behavior of the integrand one easily extracts the diver- 
gent piece of AC/4: 



A£4,<fe = -PVr+p In A [$4 + 12$ 2 J $i J 2 + 6 J $i |4j 

= S/o /3 ^W 4 (-Slnf). (A14) 



This result reproduces the usual 1-loop counterterm of the 4 -theory cor- 
rectly. 

In ( [AT3D there are 3 integrals which can be evaluated only numerically: 



h C dx[j- x m + ^ - U = & In M _ 1.4133 x 1(T 2 , 
^/o 00 ^^^T = 1-7121 xlO- 3 , 

^ Jo 00 ^iWi^T = 48875 x 10 " 4 - ( A15 ) 



The logarithmic divergence of the first integral of ( A15| ) contributes to 



( |A14| ), and is absorbed eventually by the counterterm defined through the 
renormalisation condition: 

<9 4 C/ 4| , , k , 

~^W=o = \r- (A16) 

This leads to the appearence of finite corrections to the couplings of the other 
quartic terms of the potential, which are of 0(\r) relative to their tree-level 
values: 

AC/ 4 = ^ Jd 3 x[^ (x) + 12(1 + aX^H^l 2 + 6(1 + b\ B )\^\% 

a = 3.56xl(T 3 , b = 3.17 x 1(T 3 . (A17) 

The evaluation of AUq 2 ^ is a quite tedious procedure. The easiest way to 
work out the complicated algebraic expressions of its integrand by making 
use of a symbolic programming package (MATHEMATICA). It is worthwhile 
to emphasize the importance of the correct implementation of the periodicity 
of G'( r i + r 2) on the square-interval < T\ < [3, < r 2 < (3. The result of the 
r-integrations could be expressed in analytic form, though its cumbersome 
form does not give any valuable insight. The radial k- integration has to be 
performed numerically (similarly to the case of AC/4, given in (|A13j )). The 



integrand is finite both in the infrared and the ultraviolet domain, though 
separate terms might display singularities. 
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The final expression with explicit numerical coefficients looks like 

AU 6 = §£[5.40183 i(T 6 M 6 + 2.7566 l(T 5 A 2 $ 2 |$i| 2 M 2 

+6.6586 l(T 6 A 2 |$i| 4 M 2 - 1.6991 10 5 A 3 <^>g | | 4 ] . (A18) 

This form also provides O(j^) corrections to the coefficients of lower di- 
mensional operators. All coefficients of dimension-6 operators are 0(\ 3 ). In 
addition the numerical coefficients are 2 orders of magnitude smaller than in 
front of the 0(\ 2 ) terms in AC/4. For this reason the starting values of these 
terms one can rightfully represented by (|A7|) . 

Appendix B 



1. Flow- equations for the couplings of all dimension 6 non-derivative opera- 
tors of the $ 4 -theory upon integration over $1 

The couplings appearing in the parametrisation of the action at scale K 
fllTD will flow following the equations, which are derived identically to (|2"E|). 
just including all operators up to dimension 6: 



1 


c x 2 db± 


62 x 2 




TT 2 X 2 + l + fel ' dx 


YI-k 2 x 2 +l+ 


dc 


&2C X 2 


D z x 2 


dx 


~~ 127r 2 (x 2 + l + b 1 ) 2 


IT 2 x 2 +l + bx ' 


da,2 




12D 2 X 2 


dx 


~ TT 2 {x 2 + l+b 1 ) 2 


7T 2 X 2 + l + b 1 ' 


db2 


b 2 2 X 2 


I8D4 X 2 


dx 


~ 12tt 2 (x 2 +1+6i) 2 


7T 2 X 2 + l + b! ' 


dDi 


D2C X 2 


C 3 X 2 


dx 


~ 2tt 2 (x 2 +l+b 1 ) 2 


6it 2 {x 2 +l+b 1 f i ' 



dD 2 _ (bxD2 1 cDz \ x 2 62c 2 x 2 

dx V127T 2 " 1 " TT 2 I 2 +l+bi) 2 12tt 2 (x 2 +l+6i) :i ' 

dD 3 _ ( b 2 D3 1 3£> 4 c \ x 2 b\c x ^ 

dx V 6vr 2 ' 2tt 2 > 2 +l+&i) 2 72tt 2 {x 2 +l+bi ) 3 ' 

dDi _ b 2 D 4 x 2 _ b% x 2 

dx 4tt 2 (x 2 +l+b 1 ) 2 1296tt 2 {x 2 +l+b 1 f i " 



2. Flow- equations for the couplings of all dimension 6 non- derivative 
operators of the electroeak theory upon integration over Aq 
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The extended Lagrangian density is parametrised as 



L = \F i3 F l3 + i(A0) f (A0) + §(AA)) 2 + |A W0 f + ^A 2 (K)(0t0)2 
+ \B x {K)Al + \B 2 {K){Alf + C{K)tf<j)Al + 

+ J5(J0(^V) a ^ + iW(^) + GW(ig) 3 . (B2) 

The starting (K = A) value of the couplings of the new (dimension 6) 



operators has been determined in []17] 



n _ C(3) (3a 6 3\g 4 \ 2 g 2 , 5A 3 \ F _ C(3) 9 2 / 109 9 4 , 47A<? 2 , 5A 2 a 

1024tt 4 V 16 8 4 27 ' > 1024tt 4 V 16 6 9 > 

Cfl 6 

64?r 4 ' 



F = -j&, G = 0. (B3) 



One notices that these couplings are 0(g 6 ), (here the counting A ~ g 2 is 
the convenient choice). It is assumed that this classification will not change 
to the end of the Ao-integration. Even more, we assume, that to accuracy 
0(g 6 ) G(k) = 0. All this has to be checked for selfconsistency at the end. 

With these assumptions a single Hubbard- Stratonovich transformation is 
sufficient to reach an extended Lagrangian, formally quadratic in Aq\ 

L = \F l3 F l3 + i(A0) f (A0) + §(A4>) 2 + + ^A 2 (K)(<f)^) 2 

+ WB^K) + 2^E X + 2C(K)tf</> + 2E(K)(^ ( j ) ) 2 }A 2 

+ D(K)(^^-lx 2 (B4) 

with 

a = B 2 + 2F<p ] <p. (B5) 

After integrating out the high frequency part of A Q in constant \i A)j </> 
background one finds the following gap equation for \ = 2yfax'- 

X = 2aA 2 + 3aI 1 (A), 
I n (A) = ^J k K j^, A = B x + x + 2Ctyty + 2£(0V) 2 . (B6) 

Repeating the arguments, already presented in the main text, we find the 
following contribution to the potential energy density of the reduced theory: 



AC/ = — £ ^-t (2B 2 A 2 + AFAltfcf) + 2C0V + 2E(tf ( f)) 2 ) n I n (B l ] 

4tt 2 ^ n 



(B7) 
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For the flow-equations truncated at dimension 6 operators one truncates this 
series at the third term. The following coupled set of equations is arrived at: 



dAi 3C k 2 dBi _ 3B 2 k 2 

dk TT 2 k 2 +B 1 ' dk tt 2 k 2 +B 1 

dA 2 _ 18E k 2 , 36C 2 k 2 dB 2 _ gg| k 2 

dk tt 2 k 2 +Bi ^ tt 2 (fc 2 +Bi) 2 ' dk tt 2 (k 2 +B 1 ) 2 - 

dC _ 3F k 2 , 3B 2 C k 2 

dk tt 2 k 2 +B 1 tt 2 (k 2 +B 1 ) 2 ' 

dD _ 2C 3 k 2 i 3CE k 2 

dk TT 2 (k 2 +B 1 )' i ~T TT 2 (fc 2 +Bl) 2 ' 

dE _ 6B 2 C 2 k 2 , ( 3B 2 E , 6CF \ k 2 

dk TT 2 (fc 2 +Bl)3 ^ V TT 2 ^ TT 2 > (fe 2 +Bl) 2 

d£ _ &B 2 C k 2 | 6B 2 F fc 2 

dfc TT 2 (fc 2 +Bl) 3 ^ 7T 2 (fc 2 +B!) 2 ' 

dG _ 2B| fc 2 

dfc tt 2 (fc 2 +B 1 )»- 



(B8) 



At this stage we use the knowledge of the order of magnitude of the dif- 
ferent couplings partly based on the solution of the set of equation truncated 
at dimension 4 operators (see main text), partly (for the new operators) by 
assumptions, whose consistency should be checked at the end of the calcula- 
tion: 

A l ^O(g 2 ), B l -ml~0(g 4 ), A 2 ~ 0(g 2 ), 
B 2 ~0(g A ), C~0(g 2 ), D,E,F~0(g 6 ). (B9) 

One instantly recognizes, that the right hand side of the equation for G is 
~ 0(g 12 ), therefore the assumption made for it is self-consistent. Similarly 
the right hand sides of the equations for E and F also vanish to 0(g 6 ), 
therefore they stay with their starting 0(g 6 ) values. Consistently omitting 
on the right hand sides all terms whose order of magnitude is smaller than 
0(g 6 ), one obtains the following simplified set of equations: 

dAi _ 3C k 2 dBi _ 3B 2 



dk TT 2 k 2 +B! ' dk TT 2 k 2 +B L 

dA 2 _ 36E k 2 1 36C 2 fc 2 dB 2 _ 3B 2 k 2 

dk tt 2 k 2 +B x "t" tt 2 (fc 2 +Bi) 2 ' dk tt 2 (k 2 +B 1 ) 2 - 

dC _ 3F k 2 , 3B 2 C k 2 

dk tt 2 k 2 +B! ~r ,r 2 (fc 2 +Bi) 2 ' 

dD _ 2C 3 k 2 
dk 



(BIO) 



with E, F being constants. 
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It is important to notice, that the equations for B\ and B 2 are the same 
as before, and the solution ( |30"D is accurate to 0(g 6 ). Then, one easily casts 
the equation of C(k) into the form: 



d , C , 3F k 2 



(Bll) 



dk K B 2 ! 7T 2 k 2 + Bi(fc) B 2 (jfe)' 
Its solution is simply 

r A n 2 1 
C{k)^Q 1 B 2 {k) + 6 -^B 2 {k) P (B12) 

7H 7fe p 2 + B^p) B 2 (p) 

The value of Qi is determined from the values taken at A, therefore Q% equals 
to its value determined before. 

Similarly, one transforms the equation of A\ and A 2 into 



dA 2 _ -1 o dB 2 (ri _i_ 3FB 2 rA rfpp 

rffe ~~ ±Z dk ^ TT 2 Jfe „2 + B] 
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p 2 +B 1 (p) B 2 (p)J tt 2 fc 2 +Bi' 

dM — dBifn I 3F rA rfpp 2 1 \ / R1 q\ 

dfc ~~ dfc W 1+ tt 2 Jfc p 2 +B 1 (p) B 2 (p)J- \ Dlo J 



The running of these three couplings is influenced by the dimension 6 op- 
erators only additively as it is shown by the next explicit formulae. These 
equalities are 0(g 6 ) accurate solutions of the above flow equations. They 
demonstrate explicitly the stability of the 0(g A ) solution appearing in fl52"|): 

C(0) - C(A) = Q 1 (B 2 (0) - B 2 (A)) + ^f§M J A 
A 1 (0) - A^A) = Q 1 (B 1 {0) - Bx(A)) - /o A 

A 2 (0) - A 2 (A) = 1205(^(0) - B 1 (A)) - ^ / A (B14) 

Acknowledgement 

This research has been supported by the bilateral Franco-Hungarian Bal- 
aton Research Agreement (Project No. 67). 

References 

[1] T. Appelquist and R. Pisarski, Phys. Rev. D23 (1981) 2305 
[2] N.P Landsman, Nucl. Phys. B322 (1989) 498 



29 



[3] T. Reisz, Z. f. Physik C53 (1992) 169 

[4] L. Karkkainen, P. Lacock, B. Petersson and T. Reisz, Nucl. Phys. B395 
(1993) 733 

[5] E. Braaten and A. Nieto, Effective Field Theory Approach to High Tem- 
perature Thermodynamics, NUHEP-TH-95-2, |hep-ph/9501375 



[6] for a recent review of the reduction in lattice QCD, see F. Karsch, Nucl. 
Phys. B34 (Proc. Suppl.) (1994) 63 

[7] N.F.Nicoll, T.S. Chang and H.E. Stanley, Phys. Lett. A57 (1976) 7 

[8] F.Wegner and A.Houghton, Phys. Rev. A8 (1973) 401 

[9] A. Hasenfratz and P. Hasenfratz, Nucl. Phys. B270 (1986) 685 

[10] A. Margaritis, G. Odor and A. Patkos, Z. f. Physik C39 (1988) 109 

[11] N. Tetradis and C. Wetterich, Nucl. Phys. B383 (1992) 197 

[12] S.-B. Liao and J. Polonyi, Ann. Phys. (N.Y.) 222 (1993) 122 

[13] T. Morris, Phys. Lett. B334 (1994) 335 

[14] W. Buchmiiller, Z. Fodor, T. Helbig and D. Walliser, Ann. Phys. (N.Y.) 
234 (1994) 260 

[15] J.R. Espinosa, M. Quiros and F. Zwirner, Phys. Lett. B314 (1993) 206 

[16] K. Kajantie, K. Rummukainen and M. Shaposhnikov, Nucl. Phys. B407 
(1993) 356 

[17] A. Jakovac, K. Kajantie and A. Patkos, Phys. Rev. D49 (1994) 6810 

[18] K. Farakos, K. Kajantie, K. Rummukainen and M. Shaposhnikov, Nucl. 
Phys. B425 (1994) 67 

[19] A. Ringwald and C. Wetterich, Nuvc. Phys. B334 (1990) 506 

[20] For a careful review of magnetic screening, see: A. Rebhan, Non-Abelian 
Debye screening in One-Loop Resummed Perturbation Theory, DESY- 
94-132 



30 



[21] A. Jakovac and A. Patkos, Phys. Lett. B334 (1994) 391 
[22] A. Linde, Rep. Prog. Phys. 42 (1979) 389 

[23] A.D. Linde and D.A. Kirzhnitz, Phys. Lett. 72 (1972) 471; Ann. Phys. 
(N.Y.) 101 (1976) 195 

[24] A. Jakovac, Reduction of the N-component scalar model at two loop 
level |hep-ph/9502TT3| 

[25] R.B. Pisarski, Nucl. Phys. B309 (1988) 476 



31 



